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ABSTRACT 


This thesis focuses on the modeling and simulation of the performance of M-ary 
modulation techniques. In addition to the most popular used coherent modulation 
schemes, forward error correction (FEC) coded signals, and fading channels are also 
considered. The channel is modeled for the two cases of fading and no fading, both with 
additive white Gaussian noise (AWGN). The effect of barrage noise interference is also 
considered. Binary phase-shift keying (BPSK), 16-PSK, and 16-QAM (quadrature 
amplitude modulation) communication systems are simulated and analyzed. For BPSK, 
both hard decision decoding (HDD) and soft decision decoding (SDD) are considered. 
For 16-PSK and 16-QAM, only HDD is examined. Simulation results and analytical 
results are compared. The results show that the modeling and simulation in SystemView, 


are in excellent agreement with the analytical results. 
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EXECUTIVE SUMMARY 


The purpose of this thesis is to model and simulate the most commonly used 
coherent systems, such as binary phase-shift keying (BPSK), 16-PSK, and 16-QAM 
(quadrature amplitude modulation) for different kinds of noisy and fading channels in 
order to develop a mechanism to evaluate the performance of these systems for 
conditions that do not easily yield to analysis. This is important because these coherent 
modulations are used in many communication systems including the sub-carriers in 
Orthogonal Frequency Division Multiplexed communications systems such as IEEE 
802.11la, g, and n WLANs and the Joint Tactical Radio System Wideband Networking 
Waveform. Therefore, in order to predict the performance of these systems, simulations 
typified by the simulations in this thesis are necessary. System View provides a method to 
simulate the performance of communication systems that are otherwise difficult to 
analyze. In this thesis, coherent modulation schemes with forward error correction (FEC) 
coded signals transmitted over fading channels are simulated. The channel is modeled for 
both fading and no fading with additive white Gaussian noise (AWGN). The effects of 
barrage noise interference are also considered. BPSK, 16-PSK, and 16-QAM 
communication systems are simulated and analyzed. For BPSK, both hard decision 
decoding (HDD) and soft decision decoding (SDD) are considered. For 16-PSK and 16- 
QAM, only HDD is examined. 


In order to verify the modeling and simulation procedure, the performance of 
coherent systems with no coding was examined first. Next, the performance of coherent 
systems with error correction coding was examined. Both hard and soft decision 
decoding was considered for binary systems, but only hard decision decoding was 


considered for non-binary systems. Finally, fading channels were also considered. 


Two major problems were found during the modeling and simulation. First, the 
total number of samples is restricted to 2°' (= 2,147,483,647) for all simulation loops, 


which reduces accuracy for high signal-noise-ratios (SNRs). Next, too much CPU 


XV 


resource is required to conduct each simulation. These restrictions are an artifact of using 
the student version of SystemView and presumably would not be as much of a drawback 


with the commercial version. 


The performance of the coherent systems examined, both with and without coding, 
all agreed very well with analytical results. Since the simulations developed have been 
verified for the AWGN case, they can now be modified to examine cases that are not 


possible to evaluate by analysis. 


XVI 


I. INTRODUCTION 


A. THESIS OBJECTIVE 

Wireless local area networks (WLANs) are increasingly used for both commercial 
and military applications. The objective of this thesis is to investigate the performance of 
systems that use different combinations of modulation order and code rates of forward 
error correcting codes such as those typically used in WLANs implemented with 
orthogonal frequency-division multiplexing (OFDM). This thesis focuses on modeling 
and simulation of the performance of M-ary modulation with forward error correction 
(FEC) coded signals for both channel fading and no channel fading with additive white 
Gaussian noise (AWGN). Both hard decision decoding (HDD) and soft decision 
decoding (SDD) are considered. Since the analysis with SDD is not exact when using 
non-binary modulation and a binary code, it is important to verify approximate analytical 


approaches by simulation. 


B. RELATED RESEARCH 

Wireless communications have ushered in a new era. Numerous studies have 
investigated the performance of OFDM for many different fading channels, such as 
OFDM over frequency-nonselective, fast Rayleigh and Ricean fading channels in [1] and 
[2], frequency-selective, Rayleigh fading channels in [3], frequency-selective, slow 
Ricean fading channels in [4], frequency-selective, slow Nakagami fading channels in 
[5], and frequency-selective, slow Nakagami fading channel with pulsed jamming in [6]. 
Theses papers derived analytical expressions for the different channel characteristics. 
These are just a few of many related studies. Distinct from the above works, in this thesis, 
we focus on the modeling and simulation of the performance of M-ary modulation with 
FEC coded signals over fading channels. The channel is modeled for the two cases of 
fading and no fading, both with AWGN. Binary phase-shift keying (BPSK), 16-PSK, and 
16-QAM (quadrature amplitude modulation) communication systems are simulated and 


analyzed. 


C; THESIS ORGANIZATION 

After the introduction, the thesis is organized into five remaining chapters. 
Following an introduction, Chapter II introduces the fundamental background of the 
digital communication system, including the modulation techniques. Fading channels are 
discussed in Chapter III. Forward error correction (FEC) is also introduced in this chapter. 
Next, Chapter IV explores the differences between simulation results and the analytical 
results. The use of SystemView software to simulate the system under investigation is 
also discussed in this chapter. Finally, Chapter V is a summary and recommendation for 


future work. 


D. REQUIRED SOFTWARE 
SystemView by Elanix, version 5.0 was used to build and conduct the simulations. 


MATLAB version 7.1.0 (R14) was used to generate the diagrams. 


I. DIGITAL COMMUNICATION SYSTEM 


Digital Communication has entered our daily lives in many different ways, for 
example, the radio and television sets in our living rooms and the networks around us. In 
this chapter, the fundamental background of digital communication, including various 


modulation techniques and fading channels are described. 


A. OVERVIEW 

Digital communication can be simply thought of as a series of processes from 
information input, encoding, modulating, demodulating, decoding, and information 
recovery. The basic digital communication model is described in Figure 1. 

With digital binary communication systems, one of two unique signals, say s,(t) 
or s,(t), which represent a bit “1” or a bit “0”, are transmitted each bit duration. 


Normally, the signals may be represented as a time sequences of symbols, where there 
are M distinct symbols, and each symbol represents g bits of information. Consequently, 
the relation between symbol alphabet size M , and q bits of information can be expressed 


as, 


q=log,M (2.1) 


or 


2’=M (2.2) 
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Figure 1. Digital Communication Model (After Ref. 7) 


Shannon’s theorem says that for an arbitrarily small probability of error, the 
maximum possible bandwidth efficiency is limited by the noise in the channel [8]. For a 


channel with AWGN, the capacity can be stated as [9] 


C=Wlog, [i+ 2e) (2.3) 


n 


or 





(2.4) 
where C is the channel capacity in bits/s, P. is the average received signal power, P, is 
the average received noise power, E, is the average received signal energy per bit, W is 


the noise equivalent system bandwidth in Hz, and N,/2 is the noise power spectral 


density (PSD) in W/Hz. 


For a certain E, / N,, Shannon’s theorem states that if C = R,, where R, is the bit 
rate, then it is possible to achieve error free communication. On the contrary, if R, >C, 


then error free communication is not possible. For a general expression, Equation (2.4) 


was modified as follows: 
Pe (2 1) (2.5) 


This equation is typically plotted on the bandwidth efficiency plane ( R,/W 
versus E,/N,) [11] as shown in Figure 2. If R, is the data rate in bits per second and W 


is the bandwidth occupied by the modulated radio frequency (RF) signal, then bandwidth 


efficiency 7, is expressed as [10] 


Np = (2.6) 


=| 


According to this equation above, there are two different regions on this 


bandwidth-efficiency plane. The bandwidth-limited region has 77, >1. The power-limited 


region has 7, <1. 


R/W (bit/s/Hz) 
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Figure 2. Bandwidth-efficiency plane (From Ref. 8) 
B. MODULATION TECHNIQUES 


1. 


Binary Phase-Shift Keying (BPSK) 
In binary phase-shift keying modulation, it represents a bit “1” or a bit “O” by 


changing the phase of the waveform, either zero or 1 radians (180°), respectively. The 


waveform with this sudden phase change is shown in Figure 3. 
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A general transmission signal is 


s(t) = 2A, cos[2af.t+0,(t) +8, ] (2.7) 


where 27 f, is the carrier frequency in radians/s, 0,(t) represents the data, and 0, 
is the fixed reference phase. Frequently, @,(t)= 0 is used to represent bit “1”, while 


0, (t)= 7 is used to represent bit “0”. 





Baseband Data 


waveform 
S, Sy Sy S, 


Figure 3. | Data and waveform format for BPSK (After Ref. 12) 


Zz Multi-level (M-ary) Phase and Amplitude 
a. M-ary Phase Shift Keying (MPSK) 
On the bandwidth limited region, a frequently used modulation format is 


M-ary phase-shift keying (MPSK). The transmission signal of the i" symbol is 


s(t) = 24, cos| 2.4 = ie 0 (2.8) 


where m= 1, 2,..., M. 


Unlike BPSK which only utilizes the in-phase component, MPSK uses both in- 
phase and the quadrature components. 
b. M-ary Quadrature Amplitude Modulation (MQAM) 
The constellations for both 16QAM and 16PSK are shown in Figure 4. In 
MQAM, amplitude and phase-shift keying are combined to transmit several bits per 
symbol (in this case g = 4 bits/symbol). 





Figure 4. The constellation of 1}6QAM and 16PSK (After Ref. 12) 


In this chapter, digital communication systems were introduced and modulation 
schemes such as BPSK, MPSK and MQAM were discussed. In the next chapter, the 
primary sources of channel interference are considered and error correction coding is 


introduced. 


Hil. FADING CHANNELS AND ERROR CORRECTION CODING 


Additive white Gaussian noise (AWGN) is always present in all communication 
systems and affects each transmitted symbol independently [11]. In addition to AWGN, 


the focus of this chapter is on fading channels and error correction coding. 


A. FADING CHANNELS 

Reflection, diffraction, and scattering are the main physical factors that affect the 
signals to be received. Due to these factors radio frequency channels are extremely 
random and are not easily analyzed. Many wireless communication channels do not have 
a line-of-sight (LOS) signal path. When there is no LOS, the signal is only received by 
multipath; that is, there are multiple signal paths from the transmitter to the receiver as a 
result of reflection, diffraction, and scattering of the original signal.[11] Multipath is 


illustrated in Figure 5. 





Figure 5. Multipath propagation 


Multipath can cause large fluctuations in the received signal’s amplitude, phase, 
and angle of arrival. Figure 6 represents an overview of how fading-channel can be 
categorized. Large-scale fading relates received signal strength for large changes in 
transmitter and receiver separation distance (usually much bigger than signal’s 
wavelength). Small-scale fading (changes in separation distances on the order of a half 
wavelength) is used to describe the rapid fluctuation of the amplitude of the received 


signal over a short period of time, and large-scale fading effects may be ignored. [11] 
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Figure 6. 


B. SMALL-SCALE FADING 


description 


4 
Small-scale fading due to 
small changes in position 


Time 
variance of 
the channel 
- 6 - 
Time- ; Doppler-shift 
domain Fourier domain 
Transforms description 


aw AT 18 
Z Fast Slow 
fading fading 


12 14 15 
Flat Fast Slow 
fading fading fading 


Fading-channel manifestations (From Ref. 11) 


The distance between small-scale fades is on the order of 2/2, as shown in 


Figure 7. As Figure 6 shows, the two main mechanisms that cause small-scale fading are 


time spreading of the signal due to multipath and time variance of the channel due to 


Doppler spread. 









Figure 7. 
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Small Scale 
Fading 


Large Scale 
Fading 


Distance 


Large and small scale fading (From Ref. 13) 


if Time Spreading of the Signal Due to Multipath 
As shown in Figure 6, time spreading of the signal will can be due to either 
frequency-selective fading or frequency nonselective fading (flat fading). The coherence 


bandwidth f, and the maximum excess delay T,, are two factors used to describe fading. 


m 


The reciprocal relationship between /, and 7) is given by [11] 


fo * (3.1) 


es 
sie 


where the time 7 represents the maximum excess delay, or the difference of arrival 


times of the first and the last received signal components. 


Typically, the received signal consists of several different multipath components, which 


causes the average received power S(r) to present multiple isolated peaks, as illustrated 


in Figure 8. [11] 


S(t) 


0 
meee 
Maximum excess delay 
Figure 8. Multipath intensity profile (After Ref. 11) 


The multipath intensity profile can change significantly for channels with the same value 


of T. A general parameter, the root-mean-square (rms) delay spread, was used instead of 


T,,, defined [11] 


o.=Vr°-7 (3.2) 


where r is the mean excess delay and 7 is the second moment of 7. 


There is no fixed relationship between f,, the coherence bandwidth, and o,, the 


root-mean-square (rms) delay spread. If the coherence bandwidth is defined as the 
frequency interval over which the frequency correlation function has a correlation of at 


least 0.9, the coherence bandwidth is approximately [10] 


1 
=~ —— 3.3 
to 500, (3.3) 


If the frequency correlation function is above 0.5, then the coherence bandwidth is 


empirically obtained. [14] 
IS (3.4) 


a. Frequency-Selective Fading 


For the frequency-selective case, if the coherence bandwidth f, is less 


than the symbol rate 1/7, ; that is [11] 


fo<— WwW (3.5) 


which means 
Teed, (3.6) 
where the symbol rate 1/T, is equal to the signaling rate. Figure 9(a) shows that various 


spectral components of the transmitted signal will be affected differently. 


b. Flat Fading 


Frequency-nonselective or flat-fading degradation appears whenever [11] 
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fi >W (3.7) 


or, equivalently, 


(ge (3.8) 
Figure 9(b) illustrates that all of the spectral components of the transmitted signal are 
affected in approximately the same way. 


Transmitted 


signal —___ 
W fo 
Spectral 


density 


“~ Channel 
frequency-transfer 
function 


Frequency 
(a) Typical frequency-selective fading case (fo < W) 











Channel 
Spectral frequency-transfer 
density function 
Transmitted 
signal ——> 
Ww 
Frequency 


fy | 


(b) Typical flat-fading case (fy > W) 


Figure 9. | Typical cases for frequency-selective and flat-fading (From Ref. 11) 


Zz. Time Variance of the Channel Due to Doppler Spread 


There are two parameters used to describe the time-variant nature of the channel: 
Doppler spread f, and coherence time T, . 


The magnitude of the frequency shift f, is given by [11] 


fa= 


|X 


(3.9) 
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where V is relative velocity and 4 is the signal wavelength. When the transmitter and 


receiver comes to each other, f, 1s positive, and when they move away from each other, 


f, 1s negative. 
The Doppler spread f, and the coherence time 7, are reciprocally related. [11] 


1 


T, = — 
0 ce 


(3.10) 


a. Fast Fading 
A channel is said to be fast fading if the symbol rate 1/7, (approximately 


equal to the bandwidth W ) is less than the fading rate 1/7, (approximately equal to f,); 


that is, fast fading is characterized by [11] 


W <f; (3.11) 


or 


T,>T, (3.12) 


In practice, fast fading only occurs for low data rates [11]. 


b. Slow Fading 
In the frequency domain, this implies that the Doppler spread of the 
channel is much less than the bandwidth of the baseband signals. That is [11] 


Werf, (3.13) 


or 


eyo (3.14) 
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3. Summary of Small-Scale Fading 

Depending on the relation between the signal parameters (such as bandwidth, 
symbol period, etc.) and the channel parameters (such as rms delay spread and Doppler 
spread), different transmitted signals will experience different types of fading. Signal 
time spreading (signal dispersion) and the time-variant nature of the channel are two 
mechanisms that result in small-scale fading. Figure 10 summarizes these two small-scale 


fading mechanisms. [11] 


Time-spreading Time-variant 
mechanism mechanism 
(due to multipath) (due to motion) 
FREQ-SELECTIVE FADING FAST FADING 
(Multipath Delay spread > Symbol time) (Channel fading rate > Symbol rate) 
Time-delay Dual Doppler-shift 
domain mechanisti domain 


FLAT FADING SLOW FADING 
(Multipath Delay spread < Symbol time) (Channel fading rate < Symbol rate) 
FREQ-SELECTIVE FADING FAST FADING 
(Channel coherence BW < Symbol rate) (Channel coherence time < Symbol time) 


Frequency Dual Time 
domain mechanism domain 
FLAT FADING 


SLOW FADING 
(Channel coherence BW > Symbol rate) Channel coherence time > Symbol time 


Figure 10. Summary of Small-scale fading (After Ref. 11) 


c RAYLEIGH AND RICEAN CHANNELS 

1, Rayleigh Fading Channel 

In the Rayleigh fading case, there is no LOS signal transmitted between 
transmitter and receiver and all of the received signal power is due to multipath. For 
Rayleigh fading channels (a special case of the Ricean fading channel), the received 


amplitude is modeled as a Rayleigh random variable with the probability density function 


(pdf) given by [14] 





sls Fam) (3.15) 


fs) 


where o is the rms value of the received voltage signal before envelope detection and 


o° is the time-average power of the received signal before envelope detection. 


2. Ricean Fading Channel 

When there is a dominant signal component present, such as a LOS signal, the 
effect of a LOS signal arriving with many weaker multipath signals gives rise to the 
Ricean distribution. The pdf of the Ricean distribution is given by [14] 


p(r)=—5 ox) eS) ,(A>0,r>0) (3.16) 
Oo Oo Oo 





The parameter A denotes the peak amplitude of the dominant signal and /,(¢) is the 
modified Bessel function of the first kind and zero-order. 

The Ricean distribution is often described in terms of a parameter ¢ (Ricean factor) [14], 
which is defined as the ratio between the deterministic LOS signal power and the non- 


LOS multipath signal power, given by ¢ = A’ /(20°), or in dB, 


2 





¢(dB) =10 log = dB (3.17) 


o 
As A—0,¢ —-oo(dB), and as the dominant path decreases in amplitude, the 


Ricean distribution degenerates to a Rayleigh distribution. 


D. SUMMARY OF MULTIPATH FADING CHANNELS 

The characteristic of the channel was described and multipath fading channels 
were discussed. The mechanisms and types of small scale fading channels, frequency- 
selective fading, flat fading, fast fading, and slow fading were also introduced. In this 
thesis, two widely used models for fading channels, the Rayleigh fading channel and the 


Ricean fading channel, are assumed. In the next section, error control coding is discussed. 
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E. ERROR CORRECTION CODING 

Error correction coding has made digital communications practical and popular. 
By using the error correction coding method, digital systems can easily detect and correct 
errors. All communication systems should be designed to combat the degenerative effects 


of the communication channel such as noise, interference, and fading. 


There are two basic error control strategies: automatic repeat request (ARQ) and 
forward error correction (FEC) coding. Forward error correction coding consists of 
adding a certain number of redundant bits to the actual data bits in a particular pattern 


such that recovery of the actual data bits is enhanced. [11] In this thesis, FEC is used. 


1. Forward Error Correcting (FEC) Coding 

Forward error correction (FEC) coding consists of adding a certain number of 
redundancy bits to the actual data bits in a specified pattern such that recovery of the 
actual data bits is possible. The two primary types of forward error correction codes are 


block codes and convolutional codes.[11] In this thesis, convolutional codes are used. 


In a system utilizing FEC coding, for every k information data bits, n coded bits 


are transmitted where n> k, the relation is, 


nt, =kT, (3.18) 
where 7, is the duration of a coded, or channel, bit and 7, is the duration of a data bit. 


Therefore, 


i ap oop (3.19) 
: n 


R, =— (3.20) 


where the coded, or channel, bit rate is R,. = 1/ T,. and r=k/n is the code rate. 


Since r <1, the coded bit rate is higher than the uncoded bit rate; and for a constant data 


bit rate, the addition of FEC coding increasing the bandwidth by the factor n/k =1/r. 
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[11] Also, the average energy per coded bit is less than the average energy per uncoded 
bit; and for a fixed average energy per data bit, the addition of FEC coding increases the 


probability of channel bit error. [11] 


The average transmitted power is the same whether coded or uncoded bits are 


transmitted: 


G21) 


Therefore, 


B= 78, (3.22) 


Coding gain is defined as the difference in the signal-to-noise ratios required by a 
FEC coded communication system and the same uncoded communication system to 


achieve a specific probability of bit error: 


E E 
G,=|— apes (3.23) 
N N 
0 uncoded gp 0 coded jz 


2. Hard Decision Decoding (HDD) and Soft Decision Decoding (SDD) 

In the process of decoding, two different types of “decisions” can be made in 
order to recover the original information. A “hard decision” can be made in which a 
received bit is demodulated as being either a binary one or zero. This determination is 


then sent to the Viterbi decoding algorithm for decoding of the actual information bits. 


A “soft decision” may also be made in which the decoder allows for varying 
degrees of certainty that a received bit is a one or a zero. This quantization enables the 
decoder to classify a received bit as a “strong” binary one or zero or marginally declaring 
the bit a binary one or zero. This approach improves the ability of the decoder to 


accurately estimate a received bit’s identity. [11] 


In this chapter, fading channels and error correction coding were both discussed. 
In the next chapter, the differences between simulation results and the analytical results 


are compared. 
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IV. COHERENT SYSTEM SIMULATION IN SYSTEMVIEW 


A. PERFORMANCE OF COHERENT SYSTEMS WITH NO CODING 

iL Coherent Detection (BPSK, MPSK, MQAM) 

In this subsection, the probability of bit error for binary phase-shift keying 
(BPSK), M-ary phase-shift keying (MPSK), and M-ary quadrature amplitude modulation 
(MQAM) waveforms transmitted over slow, flat fading channels is listed. The following 


equations and tables are summarized after [15]. 


For each of these signaling techniques, the conditional probability of bit error is of 


the form [15] 


Pirs) = Ol vors (4.1) 


where Q(¢) is the Q-function and a, b, and gq are constants that depend on the 
modulation type. The constants a, b, and q for BPSK, MPSK, and MQAM are listed in 
Table 1. For MQAM with qg even, a square constellation is assumed, while for MQAM 
with qg odd, a rectangular constellation is assumed. For example, in Figure 11, M =8 
(q=3, blue line) yields a rectangular constellation and M =16 (q=4, red line) has a square 


constellation. 
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Figure 11. MQAM signal constellation (After Ref. 7). 


Table 1. Modulation dependent constant for Equation (4.1) (After Ref. 15). 














Modulation q a b 
BPSK 1 1 
QPSK 2 2 1 
MPSK_ | log, M 2 2-sin?(z/M) 





MQAM_ |4,6,8,---| 4-—27%7) | 3/(2%-1) 




















MQAM_ ]3,5,7,-"° 4 3/(24-1) 





Zi Model the Coherent System (BPSK, MPSK, MQAM) in System View 
a. Simulation of BPSK in SystemView 
Baseband BPSK is the best binary modulation technique in terms of BER 
for channels with AWGN. In SystemView, a PN sequence generator is used to simulate 
the data source. Figures 12 and 13 show the simulation results for baseband BPSK, which 


shows good agreement with the analytical result in equation (4.1). 
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Performance of BPSK in AWGN 
with no channel fading. 


Final Yalue Sink 10: 

Loop N Time Value 
19080 19.0800000000000e-3 7.861635220125795e-2 
57735 57.7350000000000e-3 3.88048117966628e-2 
178432 178.432000000000e-3 1.24278151072520e-2 
806326 806.326000000000e-3 2.38893826027960e-3 
8759138 8.75913800000000 1.88612531014187e-4 
75868001 75.8680010000000 3.62098222406182e-6 
142976864 142.976864500000 1.49011614158923e-8 





Figure 12. The simulation of baseband BPSK in SystemView. 


Performance of simulated baseband BPSK in AWGN without coding 
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Figure 13. The differences of performance between simulated and theoretical BPSK. 
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b. Simulation of 16QAM in SystemView 
Figures 14 and 15 show the simulation result for 16QAM. Again, there is 
excellent agreement between the simulation results and the analytical results in equation 


(4.1). 





Performance of 16QAM in AWGN 
with no channel fading. 


Final Value Sink 8: BER 

Loop N Time Value 
10729 1.07350000000000 1.39807997017429e-1 
26470 2.64820000000000 9.52925481227368e-2 
52626 5.26440000000000 5.73482183820156e-2 
107660 10.7684000000000 2.72558781843951e-2 
276925 27.6955000000000 8.86184385431129e-3 
1104119 110.415500000000 1.81335938123342e-3 
2446290 244.633200000000 1.37091324428854e-4 
3788461 378.850927000000 2.23518463742697e-6 





Figure 14. The simulation of 16QAM in SystemView. 


Performance of simulated 16QAM in AWGN without coding 
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Figure 15. The differences of performance between simulated and theoretical 16QAM. 
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c. Simulation of 16PSK in SystemView 
Figures 16 and 17 show the simulation result for 16PSK which again 


compare very well with analytical results in equation (4.1). 





Bit error probability, P,, 


Figure 17. 


| 
0 2 4 6 8 10 12 14 16 18 


Performance of 16PSK in AWGN 
with no channel fading. 


| Final Value Sink 18: 
Time Value 
5907 591.200000000000e-3 1.69319336268202e-1 
13529 1.35390000000000 1.31216375803700e-1 
23695 2.37100000000000 9.83767830791933e-2 
38386 3.84060000000000 6.80735194009530e-2 
63947 6.39720000000000 3.91236306729265e-2 
111362 11.1392000000000 2.109081 70582528e-2 
254397 25.4432000000000 6.99134471524253e-3 
925479 92.5520000000000 1.35750921645939e-3 
1596561 159.660863000000 1.31131515969732e-4 
0 2267643 226.769726000000 2.98026172658483e-6 


Figure 16. The simulation of 16PSK in SystemView. 


Performance of Simulated 16PSK in AWGN without coding 
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The differences of performance between simulated and theoretical 16PSK. 
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B. PERFORMANCE OF COHERENT SYSTEMS WITH CODING 
1. Convolutional Codes with Hard Decision Decoding 
When a code rate r=k/n convolutional code is employed, the probability of bit 
error is upper bounded by [16] 
is <2 > BP, (4.2) 


d= free 


where d,,, is the free distance of the convolutional code, B, is the total number of 


free 
information bit ones on all weight d code sequences, where d represents the number of 


code bit one in the sequence, and P, is the probability of selecting a weight d output 


sequence as the transmitted code sequence assuming the all zeros sequence was 


transmitted. The free distance of the convolutional code d,,, is defined as the weight of 


free 


the minimum weight of non-zero code sequence. The quantities B, and d,,, are 


free 


parameters of the convolutional code chosen, and P, is determined by the type of 


modulation, the channel, and whether hard decision decoding (HDD) or soft decision 


decoding (SDD) is used. For HDD [16] 


. d 1 d-i 
& |, |ea-p) for d odd 
p — Jr (4.3) 
Saas a (d 
= [pd-p)\"" + ' Jora- p)* for d even 
2 ue a l 


where p is the probability of channel bit error. 


The free distances and the first eight B,s of commonly used r=1/2 


convolutional codes are listed in Table 2 [9], and the corresponding generator 
polynomials are listed in Table 3 [17]. The generator polynomials are a concise way of 
specifying the connections from convolutional encoder shift register input and output to 
the coded outputs. For example, the generator polynomials (131,171) (expressed in octal) 


define the r =1/2 convolutional encoder shown in Figure 18. [18] 
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Table 2. Best (maximum free distance) rate 1/2, convolutional code information weight 
structure (From Ref. 9) 
























































Vv d, ree B, free B, free tl B, free +2 B, free +3 Bu, +4 B, free +5 B, free +6 By +7 
3.5 1 4 12 32 80 192 448 1024 
4| 6 2 a 18 49 130 333 836 2069 
Dl oe 4 12 20 72 220 500 1324 | 3680 
6| 8 2 36 32 62 332 701 2342 | 5503 
7 | 10 36 0 211 0 1404 0 11633 0 

8 | 10 2 22 60 148 340 1008 | 2642 | 6748 
2 ale 35 0 281 0 2179 0 15035 0 





Table 3. Best (maximum free distance) rate 1/2, convolutional code generators (in octal) 
(From Ref. 17). 


















































v | Generators [jee] Bay | Bajees | Bajere | Bag +7 
3 7,5 5 a 192 448 1024 
4 17,15 6 2 333 836 2069 
5 3 O23 7 4 500 1324 | 3680 
6 pres: 8 2 701 2342 | 5503 
q 171,133 10 36 0 11633 0 

8 371,247 10 2 1008 | 2642 | 6748 
9 753,561 12 30 0 15035 0 
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Output 








Output 





Figure 18. The r=1/2, constraint length 7 convolutional encoder defined by the generator 
polynomials (131,171). (From Ref. 18) 


2. Simulation of Convolutional Coded BPSK with Hard Decision 
Decoding 


Figures 19 and 20 show the simulation results for BPSK in AWGN with r=1/2, 
v =3 convolutional source coding and HDD, which compare well with the analytical 


upper bound in equation (4.2) for higher SNR. 


Baseband BPSK in AWGN with y=1/2, 
¥ = 3 convolutional source coding and 
Hard Decision Decoding(HDD). 


is 


Final Value Sink 20: coded BER 


AME j|Loop N Time Value 
2313 2.34600000000000e-3 2.16169476869866e-1 
6632 6.69800000000000e-3 1.15767538782125e-1 
Sa 14655 14.7540000000000e-3 6.23208276205908e-2 
ace} 28963 29.0950000000000e-3 3.49454850433324e-2 
71811 71.9760000000000e-3 1.16691560866318e-2 
223730 223.928000000000e-3 3.29122756205610e-3 


z om” 1003754 1.00398500000000 6.41005917766633e-4 
[Ss] ef 2681442 2.68170600000000 9.95417503135267e-5 





Figure 19. The simulation of baseband BPSK in AWGN with r=1/2 , v=3 convolutional 
source coding and HDD. 
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Performance of BPSK in AWGN with r = 1/2, v = 3 convolutional source coding and Hard Decision Decoding 
T T T T T T T T 
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Figure 20. Performance for BPSK with r=1/2 , v=3 convolutional source coding and 
HDD. 


3. Simulation of Convolutional Coded BPSK with Hard Decision 
Decoding and Barrage Noise Jamming 

The primary motivation for simulating communication systems is to provide an 

evaluation tool for systems that are difficult to analyze. In this subsection, this principle is 

illustrated by using the example of BPSK with r=1/2 , v=3 convolutional coding 


where the signal is affected by barrage noise as well as AWGN. Since the barrage noise 


is just more AWGN, we expect the BER performance for cases where E,,/ N, << E,/N, 
to be approximately the sum as for the corresponding E,/ N, in the AWGN only case, as 
illustrated in Figure 20. The simulation results for E,/N,=13.35 dB and HDD are 


shown in Figures 21 and 22 and correspond well with what is expected when 


E, iN, <<£,/N,- 
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Baseband BPSK (Eb/No=13.35 dB) 
in AWGN with r=1/2, v=3 
convolutional source coding, 
bartage noise jamming, and 

Hard Decision Decoding(HDD). 


Final Yalue Sink 14: coded BER 
Loop N Time Value 
1 4804 4.83700000000000e-3 2.08159866777685e-1 
13256 13.3220000000000e-3 1.18315191670611e-1 
27462 27.5610000000000e-3 7.03927917781219e-2 
56456 56.5880000000000e-3 3.44898944609230e-2 
139333 139.498000000000e-3 1.20660738202396e-2 


1983923 1.98415400000000 6.57844965676939e-4 

3661611 3.66187500000000 9.17929913070845e-5 

5339299 5.33959650000000 1.13251093171078e-5 
7016987 7.01731800000000 1.19211677022188e-6 


2 
3 
4 
a 
3) 463808 464.006000000000e-3 3.08190153324601e-3 
rd 
8 
3 
u 
1 8694675 8.69503950000000 0.00000000000000e+0 





Figure 21. The simulation of baseband BPSK in AWGN with r=1/2 , v=3 convolutional 
source coding, HDD, and barrage noise jamming. 
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Figure 22. Performance for BPSK with r=1/2 , v=3 convolutional source coding, HDD, 
and barrage noise jamming. 
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4. Convolutional Codes with Soft Decision Decoding 
For soft decision detection with BPSK, we assume that the receiver is a maximal 
ratio combiner.[7] The probability of bit error for BPSK in AWGN with rate r =1/2 


convolutional source coding and soft decision decoding is of the form [15] 


1 | 2(1+¢) J eof -brdcy, | as 


P, = = = 
“are | 2(14+¢)+bry, 2(1+¢)+bry, 





where b is given in Table 1, c=1.0+0.1¢ is empirically obtained, ¢ is the Ricean factor, 
and y, is the average E,/N,. 


5. Simulation of Convolutional Coded BPSK with Soft Decision 
Decoding 


Figures 23 and 24 show the simulation results for BPSK in AWGN with r=1/2 , 
v =3 convolutional source coding and SDD, which compare well with the analytical 


upper bound in equation (4.2) for high SNR. For E,/N,> 6 dB, simulation results were 


not obtained due to the excessive runtime required. 


| Baseband BPSK in AWGN with r=1/2, 
¥ = 3 convolutional source coding and 
Soft Decision Decoding(SDD). 


Final Value Sink 20: coded BER 
rte) on | Time Value 
9341 9.37400000000000e-3 2.35849056603774e-1 
32519 32.5850000000000e-3 8.71080139372822e-2 
- 95520 95.6190000000000e-3 3.57568533969011e-2 
375716 375.848000000000e-3 1.25628140703518e-2 
1908536 1.90870100000000 4.17362270450751e-3 
3586224 3.58642200000000 4.86381322957198e-4 


23 7 150 5263912 5.26414350000000 6.55664223622032e-5 
es el i in 





Figure 23. The simulation of baseband BPSK in AWGN with r=1/2 , v=3 convolutional 


source coding and SDD. 
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Performance of BPSK in AWGN with r = 1/2, v = 3 convolutional source coding and Soft Decision Decoding 
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Figure 24. Performance for BPSK with r=1/2 , v =3 convolutional source coding and 
SDD. 


6. Simulation of Convolutional Coded BPSK with Soft Decision 
Decoding, and Barrage Noise Jamming 


Analogous to section IV.B.3, in this subsection the effect of both barrage noise 
and AWGN on a system utilizing BPSK with r=1/2 , v=3 convolutional coding is 
simulated, only here SDD is used. As with HDD, simulation results are obtained for 


E,/ N, =13.35 dB. We except the BER performance for cases where E,/N, << E,/N, 
to be approximately the sum as for the corresponding E, /N, in the AWGN only case, as 


illustrated in Figure 24. The results are shown in Figures 25 and 26 and, as with HDD, 


correspond well with what is expected when E,/N,>>E,/N,. 
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Baseband BPSK (Eb/No=13.35 dB) 
in AWGN with r=1/2, »=3 
convolutional source coding, barrage 
noise jamming, and Soft Decision 
Decoding(SDD). 


Final ¥alue Sink 20: coded BER 
ys) on | Time Value 
10020 10.0530000000000e-3 9.98003992015968e-2 


34698 34.7640000000000e-3 4.05219223599968e-2 
101941 102.040000000000e-3 1.48714364320450e-2 
357211 357.343000000000e-3 3.91742077016492e-3 
1874024 1.87418900000000 6.59277049972541e-4 
3551712 3.55191000000000 9.17929913070845e-5 
5229400 5.22963150000000 7.15270062133126e-6 
6907088 6.90735300000000 0.00000000000000e+0 
8584776 8.58507450000000 0.00000000000000e+0 


COnnonewn 





Figure 25. The simulation of baseband BPSK in AWGN with r=1/2 , v =3 convolutional 
source coding, SDD, and barrage noise jamming. 
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Figure 26. Performance for BPSK with r=1/2 , v=3 convolutional source coding, SDD, 
and barrage noise jamming. 
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C. FADING CHANNELS 
iL; Performance of Coherent System with Flat, Slow Fading Channels 


For BPSK received through a flat, slow, Rayleigh fading channel, the bit error 


ratio (BER) is, [15] 
a ee eee: 
P at Pe (4.5) 


For MPSK and Rayleigh fading, the BER is, [15] 


eae 
ps ee. (4.6) 
qd 1+ qy,sin’(z7/M) 


For MQAM with q even and Rayleigh fading, the BER is, [15] 





and y, is the average E,/N,. 


ps9) {| (4.7) 
q 2(2% -1)+3qy7, 


For MQAM with q odd (qg =3) and Rayleigh fading, the BER is, [15] 








P= 21 yf 3dr (4.8) 
q 2(2* -1)+ 347, 
For BPSK and Ricean fading, the BER is, [15] 
Bia Se (4.9) 
2Vac\y,+¢ +1 y,+¢+1 


where c =1.2+0.1¢ is empirically obtained. 
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For MPSK and Ricean fading, the BER is, [15] 








_qsin?(1/M)fy. 
Pia fp (4.10) 
qvac{l+¢+qsin’*(z/M)y, | 1+¢+qsin'(z/M)y, 
where c is empirically obtained. 
For MQAM with q even and Ricean fading, the BER is, [15] 
ee q/2 Ges = oR 
, 2 AUR 2792" DE +1) os __ 3067, are 
qnac|3qy7, + 2(2* -IN(¢+))] 3qg7, +2(2"-I(g +1) 
where c=1.3+0.05¢ is empirically obtained. 
For MQAM with g odd and Ricean fading, the BER is, [15] 
= q/2 = se 
qvac[3qy, + 2(2* -1I)(¢ +] 3q7, + 2(2% -I(e +1) 
where c=1.3+0.05¢ is empirically obtained. 
Zz Simulation of Coherent Systems in SystemView with Fading Channels 


a. Rayleigh Fading 


(1) BPSK. The performance of baseband BPSK transmitted over a 


slow, flat Rayleigh fading channel is shown in Figures 27 and 28. The results of the 


analytical expression in equation (4.5) are also shown in Figure 28 and shows excellent 


agreement with the simulated results. 
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Bit error probability, P, 


Figure 28. 


Figure 27. 
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EN 
2 


ete 
co) 


10 


do 


Ein 


Baseband BPSK in AWGN 
over a slow, flat Rayleigh 
fading channel 


Final Yalue Sink 10: 
Loop N Time 


Value 


23562 23.5620000000000e-3 6.36618283677107e-2 
84792 84.7920000000000e-3 2.33136462542742e-2 
284343 284.343000000000e-3 7.51687538523986e-3 
892026 892.026000000000e-3 2.46839223740009e-3 
2929494 2.92949400000000 7.60380332103714e-4 
5890612 5.89061200000000 2.47464111517227e-4 
9633910 9.63391000000000 8.10623942016546e-5 








Baseband BPSK in AWGN over a slow, flat Rayleigh fading channel. 
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The performance for simulated and approximate baseband BPSK in AWGN over 
a slow, flat Rayleigh fading channel. 
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(2) 16QAM. Figures 29 and 30 show the simulation results for 
16QAM in AWGN over a slow, flat Rayleigh fading channel. There is excellent 


agreement between the simulation results and the analytical results in equation (4.7). 





No Coding 16QAM in AWGN 
ovel a slow, flat Rayleigh 
fading channel 


Final Yalue Sink 8: BER 
ys) a | Time Value 


8116 812.200000000000e-3 9.45179584120983e-2 
19549 1.95610000000000 4.14421881475342e-2 
32502 3.25200000000000 1.40666760444507e-2 
45451 4.54750000000000 4.5836516424751 7e-3 
58277 5.83070000000000 1.52788388082506e-3 
71118 7.11540000000000 4.77144765721920e-4 
83944 8.39860000000000 1.43095635583115e-4 
97281 9.73290000000000 4.76985451943716e-5 





_ fin sia 


Figure 29. 16QAM in AWGN over a slow, flat Rayleigh fading channel. 
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Figure 30. The performance for simulated and approximate 16QAM in AWGN over a slow, 
flat Rayleigh fading channel. 
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(3) 16PSK. Figures 31 and 32 show the simulation results for 16PSK 


in AWGN over a slow, flat Rayleigh fading channel. There is good agreement between 


the simulation results and the analytical results in equation (4.6). 


7 13 ; . 
: 2 No Coding 16PSK in AWGN 
over a slow, flat Rayleigh 
fading channel 


18 
‘Esl 


Final Yalue Sink 18: 
bys) sa | Time Value 
721 72.6000000000000e-3 1.07642626480086e-1 


1194 120.400000000000e-3 5.50357732526142e-2 
5873 588.800000000000e-3 2.28258388495777e-2 
14719 1.47390000000000 7. 72974520469510e-3 
45434 4.54590000000000 2.48115278175398e-3 
254293 25.4323000000000 8.58860578299456e-4 
589831 58.9867000000000 2.86286859433152e-4 
925369 92.5411310000000 9.54289531443840e-5 





Figure 31. 16PSK in AWGN over a slow, flat Rayleigh fading channel. 
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Figure 32. 
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The performance for simulated and approximate 16PSK in AWGN over a slow, 
flat Rayleigh fading channel. 
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b. Ricean Fading, ¢ =4 


(1) BPSK. Figures 33 and 34 show the simulation results for 
baseband BPSK in AWGN over a slow, flat Ricean fading channel with ¢ =4. There is 


good agreement between the simulation results and the analytical results in equation (4.9). 





Baseband BPSK in AWGN 
over a slow, flat Ricean 
fading channel, ¢ = 4. 


Final Value Sink 10: 

Loop N Time Value 
32153 32.1530000000000e-3 3.11012969240817e-2 
255680 255.680000000000e-3 5.55308751665926e-3 
1562441 1.56244100000000 1.02998397802701e-3 
4917884 4.91788400000000 2.42292895453745e-4 
8273327 8.27332700000000 7.15256419426364e-5 
11628770 11.6287700000000 2.14576925827909e-5 





Figure 33. Baseband BPSK in AWGN over a slow, flat Ricean fading channel with ¢ =4. 
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Performance of baseband BPSK in AWGN over a slow, flat Ricean fading channel with ¢ = 4 
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Figure 34. The performance for simulated and approximate Baseband BPSK in AWGN over 
a slow, flat Ricean fading channel with ¢ =4. 


(2) 16QAM. Figures 35 and 36 show the simulation results for 
16QAM in AWGN over a slow, flat Ricean fading channel with ¢ =4. There is excellent 


agreement between the simulation results and the analytical results in equation (4.11). 


> — Ein” No Coding 16QAM in AWGN 
over a slow, flat Ricean 
. fading channel, € =4. 


meee CE] 


Final Yalue Sink 8: BER 

sys) sm | Time Value 
11514 1.15200000000000 6.47149460708783e-2 
32910 3.29220000000000 1.69491525423729e-2 
59491 5.95090000000000 3.06748466257669e-3 
$9357 8.93810000000000 5.72573718866304e-4 
119723 11.9753000000000 1.43095635583115e-4 
149013 14.9049000000000 4.76905834942891e-5 





4 ~ fin al 


Figure 35. 16QAM in AWGN over a slow, flat Ricean fading channel with 6 =4. 
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Performance of baseband 16QAM in AWGN over a slow, flat Ricean fading channel with ¢ = 4 
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Figure 36. The performance for simulated and approximate 16QAM in AWGN over a slow, 
flat Ricean fading channel with ¢ =4. 


(3) 16PSK. Figures 37 and 38 show the simulation results for 16PSK 
in AWGN over a slow, flat Ricean fading channel with ¢ =4. There is very good 


agreement between the simulation results and the analytical results in equation (4.10). 


— [3 No Coding 16PSK in AWGN 
Oo over a slow, flat Ricean 

fading channel, € = 4. 

—i> 


Final Value Sink 18: 

rts) a | Time Value 
10171 1.01760000000000 7.94912559618442e-2 
42470 4.24800000000000 3.69822485207101e-2 


210236 21.0252000000000 8.79541108986616e-3 
378002 37.8024150000000 1.52963671128107e-3 
545768 54.5796300000000 2.92073483304126e-4 
713534 71.3568450000000 7.74888833255844e-5 





Figure 37. 16PSK in AWGN over a slow, flat Ricean fading channel with ¢ =4. 


39 


Performance of baseband 16PSK in AWGN over a slow, flat Ricean fading channel with ¢ = 4 
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Figure 38. The performance for simulated and approximate 16PSK in AWGN over a slow, 
flat Ricean fading channel with ¢ =4. 


i Ricean Fading, ¢ =10 
(1) BPSK. Figures 39 and 40 display the simulation results for 
baseband BPSK in AWGN over a slow, flat Ricean fading channel with ¢ =10. There is 


good agreement between the simulation results and the analytical results. 
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Baseband BPSK in AWGN 
over a slow, flat Ricean 


fading channel, ¢ = 10. 


Final Yalue Sink 10: 

Loop N Time Value 
49593 49.5930000000000e-3 1.61313088540721e-2 
904344 904.344000000000e-3 8.84619127234769e-4 
26453434 26.4534340000000 3.02418204003306e-5 
26843545 26.8435450000000 2 309680036671 76e-6 
81629741 81.6297418000000 4.09781941990151e-7 
108473286 108.473287200000 1.11758711451859e-7 





Figure 39. Baseband BPSK in AWGN over a slow, flat Ricean fading channel with ¢ =10. 
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Figure 40. The performance for simulated and approximate Baseband BPSK in AWGN over 
a slow, flat Ricean fading channel with ¢ =10. 
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(2) 16QAM. Figures 41 and 42 display the simulation results for 
16QAM in AWGN over a slow, flat Ricean fading channel with ¢ =10. There is good 


agreement between the simulation results and the analytical results. 


% No Coding 16QAM in AWGN 
|@| over a slow, flat Ricean 
fading channel with € = 10. 


Final Yalue Sink 8: BER 
Loop N Time Value 
15694 1.57000000000000 4.68164794007491e-2 


60043 6.00550000000000 7.64818355640535e-3 

133773 13.3791000000000 3.33889816360601e-4 
2226339 22.2663000000000 1.19217930376725e-5 
315124 31.5154000000000 1.25607100988109e-6 
406005 40.6041000000000 2.42292895453745e-7 





Figure 41. 16QAM in AWGN over a slow, flat Ricean fading channel with ¢ =10. 


Performance of 16QAM in AWGN over a slow, flat Ricean fading channel with ¢ = 10 


T T 








ul 


I 
Vv Approximate 
—4— Simulate 























ial c a 
3 [ * ] 
= | Bi 4 
alae, 
2 10 E aN | 
a: K 1 
S \ 1 
5 2 | 
eS 40eb * J 
~ ] 
10° ee 4 
a 
10°- ' 4 
5 10 15 20 25 30 
E, /N, (dB) 


Figure 42. The performance for simulated and approximate 16QAM in AWGN over a slow, 
flat Ricean fading channel with ¢ =10. 
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(3) 16PSK. Figures 43 and 44 display the simulation results for 
16PSK in AWGN over a slow, flat Ricean fading channel with ¢ =10. There is a good 


agreement between the simulation results and the analytical results. 





‘ No Coding 16PSK in AWGN 
over a slow, flat Ricean 
fading channel, € = 10. 
[23] 


Final Value Sink 18: 

Loop N Time Value 
607 61.2000000000000e-3 6.47149460708783e-2 
2462 247.200000000000e-3 2.69687162891046e-2 
31412 3.14270000000000 3.14915545376467e-3 
1373583 137.360400000000 1.43095635583115e-4 
2715754 271.578127000000 5.50357732526142e-6 
4057925 405.795854000000 5.96057674540589e-7 





Figure 43. 16PSK in AWGN over a slow, flat Ricean fading channel with ¢ =10. 
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Figure 44. The performance for simulated and approximate 16PSK in AWGN over a slow, 
flat Ricean fading channel with ¢ =10. 
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C. CHAPTER SUMMARY 

In this chapter, the bit error ratio performance of various types of modulated 
signals, with and without error correction coding and with and without channel fading, 
were obtained by simulation and the results compared with those obtained analytically 
when possible. In the next chapter, conclusions are drawn and suggestions for future 


work made. 
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V. CONCLUSION AND FUTURE WORK 


A. CONCLUSION 

Using SystemView, we simulated the effect of various types of fading channels. 
Two major problems were found during the modeling and simulation. First, the total 
number of computations is restricted on each simulation loop, which limits the degree of 
accuracy for high SNR. Next, too much CPU resource is required for each simulation, 
and it is time consuming to accurately simulate performance for high SNR. Presumably, 
using the commercial version instead of the student version would alleviate these 


difficulties. 


Except for these two problems, generally speaking, SystemView provides a good 
modeling and simulation tool for communication systems. The performance of coherent 
systems with no coding, the performance of coherent systems with coding, and the 
performance of coherent systems with no coding and fading channels all agree very well 
with analytical results. Most importantly, the performance of communication systems 
affected by narrowband noise and other channel impairments that are difficult to analyze 
are easily simulated in SystemView. Since the simulations developed have been verified 
for the AWGN case, they can now be modified to examine cases that are not possible to 


evaluate by analysis. 


B. FUTURE WORK 

With rapid progress in hardware and software, the simulation process will be 
significantly improved. The utility of this simulation tool for evaluating difficult to 
analyze systems has been demonstrated, and SystemView can be used to evaluate any 
number of modern communication systems. The combination of error correction coding 
and modulation specified by the various digital standards, such as JEEE 802.11g, should 
be examined for fading channels and interference conditions other than simple AWGN to 
foster improved understanding of how these popular communications wireless standards 
perform in various radio environments. For example, the barrage noise jamming 


considered in Chapter IV is simply additional AWGN. It would be valuable to re- 
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examine this case, but for a barrage noise jammer implemented by frequency modulating 
a voltage controlled power oscillator using bandlimited AWGN or bandlimited AWGN 


plus a periodic waveform, as is typical of real military electronic warfare systems [19]. 
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